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Abstract. We propose integrable discretizations of derivative nonlinear Schrodinger 
(DNLS) equations such as the Kaup-Newell equation, the Chen-Lee-Liu equation 
and the Gerdjikov-Ivanov equation by constructing Lax pairs. The discrete DNLS 
systems admit the reduction of complex conjugation between two dependent variables 
and possess bi-Hamiltonian structure. Through transformations of variables and 
reductions, we obtain novel integrable discretizations of the nonlinear Schrodinger 
(NLS), modified KdV (mKdV), mixed NLS, matrix NLS, matrix KdV, matrix mKdV, 
coupled NLS, coupled Hirota, coupled Sasa-Satsuma and Burgers equations. We also 
discuss integrable discretizations of the sine-Gordon equation, the massive Thirring 
model and their generalizations. 
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1. Introduction 



The inverse scattering method (ISM) was invented by Gardner et al. [jl| more than 
thirty years ago. They expressed the KdV equation as the compatibility condition of 
an eigenvalue problem and time evolution of the eigenf unction, and solved the KdV 
equation through the inverse problem of scattering. A pair of operators which defines 
the eigenvalue problem and the time evolution is now called the Lax pair. Zakharov 
and Shabat [|2|, |^ considered a generalization of the eigenvalue problem and solved the 
nonlinear Schrodinger (NLS) equation 



iut + Uxx — ^uvu = O 



+ 2vuv = O 



'1.11 



under some conditions via the ISM. It is noteworthy that the NLS equation (|1.1| ) is 
integrable for matrix- valued variables u and v ||, ^ |^. Throughout this paper, we 
use the symbol italic O when dependent variables in the equation considered can take 
their values in matrices. Ablowitz et al. ^ formulated the Zakharov-Shabat method in 
a plain manner and constructed the hierarchy of the NLS equation (|1 . 1|) . Up to now, a 
variety of nonlinear evolution equations in continuous space-time have been solved via 
the ISM based on modified versions of the Zakharov-Shabat eigenvalue problem (see, 
e.g. |, |10|, |n|, 0, 0, |1|, [1|, |11 



In the remarkable paper |T^, Ablowitz and Ladik proposed a natural discretization 
of the Zakharov-Shabat eigenvalue problem and obtained a discrete NLS hierarchy. 
Their success suggests that discrete integrable hierarchies are obtained through natural 
discretizations of the modified Zakharov-Shabat eigenvalue problems. It is, however, 
quite difficult to find such discretizations directly. In a few successful studies [|^, pO|, [2l[], 
the difficulty is avoided skillfully with the help of alternative approaches. Hence, wide 
applicability of Lax-pair formulations has not been established in the theory of discrete 
integrable systems. 

In this paper we discuss integrable discretizations of derivative NLS (DNLS) 
equations and related systems by elaborating a new formulation of Lax pairs. Among 



I.e. 



the DNLS equations, the following three representatives are well-known [p2| , |23| , p4 
the Kaup-Newell equation [p!5[] , 

iqf + g^,. - i(gV-g«). = O irf - r^, - i(r-g«r-),. = O (1.2) 

the Chen-Lee-Liu equation p5| , 

O irf-r-^-ir-q-r- = (1.3) 



^1? + Ixx 



and the Gerdjikov-Ivanov equation p6[]. 



^Qt + Q.XX + n'^r'^q'' + ^q''r°q''r°q'' 



ir. 



_j_ * G^G^G 
XX ' 



_^G^G^G^G^G 



o 
o. 



;i-4) 



It has been shown that ( |1.2| )-( pr^ ) are integrable for matrix-valued variables |2^, ^ 
(see [^, 1^ for complete hsts of integrable coupled DNLS equations). The DNLS 
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equations (p..2|)-(pr4D correspond to special cases (5 = —1/2, —1/4, 0, respectively) of 
the generalized DNLS equation |^ 

+ Qxx + i(45 + l)gV, + AiSqq^r + {5 + 1/2) (45 + 1)^^^ = 

in - '^xx + i(45 + ly^qx + 4i5rr^g - (5 + 1/2) (45 + l)r^q^ = 0. 

The zero symbol indicates that the dependent variables are restricted to scalars. We 
generate ( |1.5| ) from the Gerdjikov-Ivanov equation ( |1.4| ) via the transformation 

g = exp (^-2i5 J g°r°da;'j r = r° exp (^2i5 J q'^r'^ dx'Y (1.6) 

A comprehensive description of physical applications of the DNLS equations as well 
as their exact solutions can be found in p3| , p^ . An integrable discretization of the 
Chen-Lee-Liu equation ( p..3|) was proposed by Date, Jimbo and Miwa |Q. However, 
their scheme does not admit the reduction of complex conjugation between and r^. 
Therefore it is of little practical use. In other studies on discrete DNLS equations 
P^ , P^l , equations of motion are not explicitly given. Practically, little is known 
about natural integrable discretizations of the DNLS equations. 

To solve this problem, we consider discrete analogues of relations between the 
DNLS hierarchies and a generalization of the NLS hierarchy |T3|, P^, HTl, 133]. We 



extend the Lax-pair formulation of Ablowitz and Ladik to obtain a generalization of 
the discrete NLS (Ablowitz-Ladik) system. We find out that integrable discretizations 
of the DNLS equations ( |1.2|) -( |0|) together with their higher symmetries are embedded 
in the generalized Ablowitz-Ladik system. With appropriate gauge transformations, we 
obtain standard forms of eigenvalue problems for the discrete DNLS systems. Using a 
discrete version of the transformation (|1.6|) , we obtain a discrete generalized DNLS 
system. All the discrete DNLS systems but one exception admit the reduction of 
complex conjugation. Thus they inherit the crucial property of the continuous DNLS 
systems. Through simple transformations and reductions for the discrete DNLS systems, 
we obtain integrable discretizations of various systems, e.g. the mixed NLS, matrix NLS, 
matrix KdV, matrix modified KdV (matrix mKdV), coupled Sasa-Satsuma and Burgers 
equations. In this way, we can derive integrable discretizations of surprisingly many 
systems from one generalization of the Ablowitz-Ladik formulation. This exemplifies 
that some approaches based on a Lax-pair formulation are very fruitful in the study of 
discrete integrable systems. 

The sine-Gordon equation and its generalization are integrable via the ISM based on 



the Zakharov-Shabat eigenvalue problem |3^, 0. In contrast, Ablowitz and Ladik were 



unsuccessful in obtaining a local discretization of the sine-Gordon equation (cf (2.10) in 
||T8|). This does not indicate that their discretization of the Zakharov-Shabat eigenvalue 
problem is useless in studying discrete sine-Gordon equations. It is just because time 



evolution of the eigenfunction assumed in |18| is too general and inappropriate. We 



prove that discrete sine-Gordon equations and their generalizations are associated with 
the Ablowitz-Ladik eigenvalue problem (see Appendix A and Appendix B ). The massive 



Thirring model and its variant types are known to have the same eigenvalue problems 
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as the DNLS equations have |T^, |3^. In particular, the massive Thirring 

model and the Chen-Lee-Liu equation ( |1.3| ) have the eigenvalue problem in common. 
This indicates that the massive Thirring-type models are related to some generalization 
of the sine-Gordon equation. Combining these items of information, we derive integrable 
discretizations of the massive Thirring-type models from the eigenvalue problems for the 
discrete DNLS systems. 

The paper consists of the following. In section 2, we show that the DNLS equations 
are embedded in a generalization of the NLS equation, and propose a generalization of 
the Ablowitz-Ladik system. In section 3, we reveal that space discretizations (for short, 
semi-discretizations) of the DNLS systems are embedded in the generalized Ablowitz- 
Ladik system. In section 4, we obtain a variety of integrable lattice systems related 
to the semi-discrete DNLS systems. In section 5, we investigate semi-discretizations of 
the massive Thirring-type models. The last section, section 6, is devoted to concluding 
remarks. 

In the present version of this paper, we do not discuss time discretizations of the 
semi-discrete DNLS systems. Some of the results on this problem can be found in old 



versions of this paper ( |arXiv:nlin. SI/010505^ ver. 1 or 2). 



2. Preliminaries 



In this section, we first show that the DNLS equations (|1.2| )-( p~^ ) are embedded 
in a generalization of the NLS equation ( |1.1| ). Next we extend the Ablowitz-Ladik 
formulation to obtain a semi-discretization of the generalized NLS system. 



2.1. Embeddings of the DNLS equations into a generalized NLS equation 

Let us begin with a brief description of Lax pairs. We consider a pair of linear equations 

= Um = V"* (2.1) 

for a column vector ^ . The subscripts x and t denote the partial differentiation by x 
and t respectively. U and V are square matrices which depend on a parameter. The 
compatibility condition of ( p.l| ) is given by 

Ut -V^ + UV -VU = (2.2) 

which we call the zero-curvature condition. If we choose the matrices U, V appropriately, 
( p. 21) gives a system independent of the parameter. In such cases, the pair of matrices U, 
V and the parameter are called the Lax pair and the spectral parameter, respectively. 
We introduce the following form of the Lax pair: 



U = iX 



V = iX' 



' -h 




w u 


+ 




h 


V s 



-2/i 



2/o 



+ A 


2u 




2v 


+ i 





— Vx — VW + SV 



wu 

d 



us 



(2.3a) 



(2.3b) 
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Here A is the spectral parameter. U and V are divided into four blocks as (/1+/2) x (^1+^2) 
matrices. Ji and I2 are, respectively, the li x li and I2 x I2 unit matrices, w, s, u, v 
are matrix-valued variables, c and d are arbitrary functions at this stage. Substituting 
( p.3| ) into the zero-curvature condition ( p.2| ), we obtain the following system (see (26) 
in for the case of scalar variables): 

iwt + Cx — wc + cw + {uv)x + uvw — wuv = O 

ist + dx — sd + ds — {vu)x — vus + svu = O {2 A) 

iut + Uxx — ud + cu — WxU — 2wUx + 2uxS + uSx + w'^u — 2wus + us"^ = O 
ivt — Vxx — vc + dv + SxV + 2sVx — 2vxW — vWx — s'^v + 2svw — vw'^ = O. 



We call (|2^ ) the generalized NLS equation. If we set w = O, s = O, c = —uv, d = vu, 
( p.4|) is reduced to the matrix NLS equation ( |1.1|) . 

We already know Lax pairs for the matrix DNLS equations ( |L2|) -( pr4D |29]. With 
the help of gauge transformations, we can transform the Lax pairs into the form (|2.3| ) 
(see, e.g. [0). Thus the DNLS equations are embedded in the generalized NLS equation 
( p.4|) . We omit details and show the main results in the following. For simplicity of the 
embedding formulae, we change the scalings of variables in ( p..2|) -( p~4D . 

(a) If q^ and satisfy the Kaup-Newell equation, 

+ ?L + 2(gV-g-). = O irf - r^, + 2(r-gV-). = O (2.5) 
w, s, u, V defined by 

yj = -q^r'' s = r''q'' u = q"" v = - r'^q'^r'' 
satisfy the generalized NLS equation (^) with 

c = -q^r"" - 2g'^r«g'^r^ d = r^'q^ + 2r«g'^r«g'^. 

(b) If q'^ and r*^ satisfy the Chen-Lee-Liu equation, 

+ f.. + '2q'.r-q- = O irf - rf^ + 2r-q-r- = O (2.6) 
w, s, u, V defined by 

w = O s = r'^q'^ u = (f = rf. 
satisfy the generalized NLS equation (|2.4| ) with 

c = -q°r1 d = r^q^ + r^q'^r^q'^. 

(c) If q° and r° satisfy the Gerdjikov-Ivanov equation, 

igf + - 2q^r^q^ - 2q^r^q^r^q^ = O 

irf - rf^ - 2r°q'^r'^ + 2r^q'^r^q°r^ = O 
u and V defined by 

u = q f = + r q r 
satisfy the matrix NLS equation ( |1.1| ). 



(2.7) 
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2.2. Generalization of the Ablowitz-Ladik system 

We consider a semi-discrete version of the linear problem ( pTl 



n,t 



The compatibility condition of (p.8|) is given by 

L„t + LnMn - Mn+lLn = O 



:2.8) 



(2.9) 



which is a semi-discrete version of the zero-curvature condition (|2.2|) . 

We generalize the Lax pair proposed by Ablowitz and Ladik [|18| and consider the 
following form: 



ZWr. 



Ur 



z'^ali + c„ 



.(2.10) 



Here z is the spectral parameter and a, 6 are constants. L„ and M„ are divided into 
four blocks as (Zi + I2) x (/i + ^2) matrices. Ji and J2 are unit matrices. w„, m„, f„ 
are matrix-valued variables. Cn and (i„ are arbitrary functions. Substitution of the Lax 
pair ( ^lOD into the zero- curvature condition ( p. 91) yields the following system: 

Wn,t + WnCn - Cn+lWn + aUnVn-lW^^-i - aW^^l^Un+lVn = O 

Sn,t + Sndn — C^n+l-^n + bVnUn-lS^^i — bs^^^^Vn+lUn = O /o i i \ 

Un,t + Undn - Cn+lUn + bWnUn-lS.^_i - aW.^_^_iUn+lSn = O 
Vn,t + VnCn - dn+lVn + aSnVn-lWn\ " bs:^liVn+lWn = O. 

We call ( |2.11| ) the generalized Ablowitz-Ladik system. When just two of w„, s„, m„, Vn 
are functionally independent, we may determine c„, dn from the consistency of ( |2.11|) . If 
we set Wn = h, Sn = h, Cn = -a{Ii + ■u„f„_i), dn = -b{l2 + VnUn-i), ( |2.11| ) is reduced 
to the matrix Ablowitz-Ladik system [|3|, |4|, g^l^ 

Un,t - aUn+l + bUn^l + {a - b)Un + aUn+lVnUn - bUnVnUn-l = O 
Vn,t - bVn+1 + aVn^l + {b - a)Vn + bVn+lUnVn " af„M„t;„_i = O. 

The system ( |2.12| ) with b = a* admits the reduction of complex conjugation Vn = crUn 
and not the reduction of Hermitian conjugation t>„ = aul^. Throughout this paper we 
use a to denote an arbitrary (but usually nonzero) real constant. We can construct an 
infinite set of conservation laws for ( p.ll ) by a recursive method [0, The first five 
conserved densities are 

trlogW^ trlogS^ tT\og{Wn- UnSn^Vn) tl{UnVn~lWn^^Wn^) tl {VnUn^lS n\Sn^) . 

The generalized Ablowitz-Ladik system (|2.11| ) gives a semi-discretization of the 
generalized NLS equation (|2.4|) together with its higher symmetry. To see this, we 
suppose the following scalings for (|2.11|) : 

Wnit) = h + Aw{x,t) Snit) = I2 + As{x,t) Un{t) = Au{x , t) Vn{t) = Av{x,t) 

a = -b=-^ Cn{t) = -^Ii + ic{x,t) dn(t) = ^l2 + id{x,t). 



(2.12) 



A2 



A2 
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Here A denotes the lattice spacing and x = nA. In the continuum hmit A ^ 0, ( |2.11| ) 



is precisely reduced to ( |2.4| ). A correspondence between the spectral parameters is given 
hj z = exp(— iAA). In the next section, we show that semi-discrete DNLS systems are 
embedded in through lattice analogues of the formulae (a)-(c) in section We 

also mention that a system of semi-discrete coupled NLS equations is obtained as a 
reduction of ( CT) g§. 

3. Semi-discrete DNLS systems 

In this section, we obtain semi-discrete DNLS systems by considering lattice analogues of 
(a)-(c) for the generalized Ablowitz-Ladik system ( |2.11|) . Through a discrete analogue 
of the transformation (|1.6|), we obtain a semi-discrete generalized DNLS system. 



3.1. Semi-discrete Kaup-Newell system 

We find that a lattice analogue of (a) is given as follows. 
(A) If and r^ satisfy the system 

+ b{h + q^^,r^)-X-i = (3.1a) 
rlt - Kh + rl^.qlVr-l^, + Kh + r^g^J-V^ - a{h - <q^)-'r^ 

+ a{h - rUqn-irX-i = O (3.1b) 
Wn, Sn, Un, f„ defined by 

Wn = /l — qnl^n Sn = I2 + '"n+l^'n '^n = qn '^n = 1^n+l ~ "^n ^ '''n+lQ'n'^n (3-2) 

satisfy the generalized Ablowitz-Ladik system (|2.11| ) with 

Cn = -hh - a{h - q^r^r' + Kh + g^iO"^ 
rf„ = -ah + a{h - <q^)-' - Kh + rlq^-i)-' ■ 
The system ( p.l| ) with a = —h = i gives an integrable semi-discretization of the Kaup- 
Newell equation ( p. 51 ), while ( p.l| ) with a = h = —1 gives an integrable lattice version of 
a higher symmetry of ( |2.5| ) (cf (3.39) in [^). 

To consider the reduction of complex conjugation or Hermitian conjugation, we 
suppose that q^ and are defined on the fractional lattice n G Z/2. If 6 = a*, the 
semi-discrete Kaup-Newell system ( p.l|) admits either the reduction = 'i<yq^ \ or the 

reduction = 'i(jq^\. This property is crucial for various applications of ( p.l|) . 
Through the gauge transformation 

f{z)h o 



z n ^2. 



(3.4) 



where ji^z) is an arbitrary function of 2, the linear problem ( p.8| ) and the Lax pair ( |2.1(]| ) 
with ( |3.2|) and (^]3|) are changed into 



$ _Li = $ * = M"^^ 
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zh 
O 



h 



o 



Ml 



[a{z^-l) + h{-l + j,)]h 
+ h{l-^){h + q^_,rl)-^ 



+ 6(-i + ^)(/2 + rX-i)-V^] 



f{z)[az{h-qlrl)-^ql 
+ f(/i + g;^_ir-)-^g-_i] 



(3.5a) 
(3.5b) 



We can check that substitution of the Lax pair (|3.5|) into the zero-curvature condition 



( p.9| ) yields the semi-discrete Kaup-Newell system ( p.l|) . If we set 

/(^) = (-l + l)^ z = exp(-iC'A) 



iA 1 



,_iA 1 

. 2 



the L„-matrix (|3.5a|) has the following asymptotic form: 

'2 ; 



' h 


+ A 


h 





Consequently, in the continuum limit A 



+ 0(A2). 

0, we recover the eigenvalue problem 



proposed by Kaup and Newell 



For the time being, we assume without loss of generality that q^ and are square 
matrices. If necessary, we append dummy variables to and which are finally 
equalized to zero. Then a set of Hamiltonian and Poisson brackets for the semi-discrete 
Kaup-Newell system (|3.1|) is given by 

H = ^[-alogdet(J - gX) + 61ogdet(J + q^r^^,)] (3.6a) 



and 



(3.6b) 



K 

qn,t 



Here 11 denotes the permutation matrix: II^-'; = 6i^i6j^k- It is easily proven that 
and r^j. = {r^,H} coincide with ( p.la|) and (|3.1b|) , respectively. 
Let us prove that in the periodic boundary case {qn+M — qni ''^n+M — ''^n) (|3-1| ) 
has an involutive set of conserved quantities. We introduce new variables by 
qn — Pn+i ~ Pn Pn+M = Pny which Satisfy ultra-local Poisson brackets: 

{Pn^pfn} = K®r-} = o K®C} = 5„,^n. 

For convenience, we set f{z) = 1/z. Applying a gauge transformation to the L„-matrix 
( p.5a| ), we obtain a new ultra-local L„-matrix: 



O I 



-I Pi 

O I 



zl 



pK 



-^nPn 
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The matrix C^{z) satisfies the fundamental r-matrix relation: 

{C(A)®C(^)} = 5„,^[C(A) ® r(A,/i)]. 
Here r(A, ii) in the case of scalar variables is given by 

1 

(A2-l)(/i2-l) ' 



(3.7) 



r(A,/i) 



A2 



1 

1 



(3i 



r(A,/i) for matrix-valued variables is given by replacing the ones in the matrix (|3.8| ) 
with the permutation matrix 11. This is similar to the r-matrix given in 0]. 

The fundamental r-matrix relation ( |3.7| ) results in commutativity of the traces of 
the monodromy matrix C\.^C\j_i ■ ■ ■ JO^ for different values of the spectral parameter. 
This generates a set of M conserved quantities in involution. Owing to the relation 
tr ■ ■ ■ £5^) = ti {Ul,jU^_^- ■ ■ L^) we can write the conserved quantities in 

terms of the original variables and r^. The set of conserved quantities includes 
the Hamiltonian (|3.6a|) itself. 



3.2. Semi- discrete Chen-Lee-Liu system 

We find that a lattice analogue of (b) is given as follows. 
(B) If and r^ satisfy the system 

Vn defined by 

h Sn = I2 + rn+lQn 



O 
O 



(3.9) 



Ur. 



n+1 



(3.10) 
(3.11) 



satisfy the generalized Ablowitz-Ladik system (|2.11| ) with 

Cn = -ah - ag°(r^ - r^_J rf„ = -6(/2 + r^qn-i]~^ 

The system ( |3.9| ) with a = —b = i gives an integrable semi-discretization of the Chen- 
Lee-Liu equation (|2.6| ), while ( |3.9| ) with a = b = —1 gives an integrable lattice version 
of a higher symmetry of ([2l^ ) (cf (3.37) in p7| ). 

Unfortunately, the semi-discrete Chen-Lee-Liu system ( |3.9| ) admits neither the 
reduction r° = icrq'^*;. nor the reduction r^ = iag^J^^. This is an intrinsic drawback for 
practical use. We obtain an alternative semi-discretization of the Chen-Lee-Liu system, 
which admits the reduction of complex conjugation, in section ^.5| . 

Through the gauge transformation ( |3.4| ) with replacing r^ by r^, the linear problem 
( p.8|) and the Lax pair ( p. 101) with (|3.10|) and ( p.llj) are changed into 

and 

zh + Iq^r^ f{z)q^ 



(3.12a) 
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[a(^2_2) + A]/^ + a(/i + gX-i) 



c^-l 

n) 



1 



nJ 



/(z)[a2;gc 



6(-l + ^)(/2 + <gS-,)-i 



(3.12b) 



We can prove that substitution of the Lax pair ( p.l2|) into the zero-curvature condition 
( ^I9D yields the semi-discrete Chen-Lee-Liu system ( p.9|) . With appropriate scahngs, 
we reproduce the eigenvalue problem for the Chen-Lee-Liu system in the continuum 
limit. 



3.3. Semi- discrete Gerdjikov-Ivanov system 

We find that a lattice analogue of (c) is given as follows. 

(C) If g° and r° satisfy the system 



+ 0'Qn+l''"n+l1n''"nQn 



G^G „G^G„G 



- H 



n \ n-\-l 



<,t - K+1 + + a)r^ - 6r^+i(g^ - q^_,)r^ + ar^{q^ - q^_,X 



Un and Vn defined by 



^^nQn^nln-l'^n-l 



o 



Ur 



/V.G1 G I G n^r^ 

' n+1 ' n ' ' n+l'in ' n 



(3.13a) 
(3.13b) 

(3.14) 



satisfy the matrix Ablowitz-Ladik system ( p.l2|) . 



If g° = /, the relation ( p.l4| ) between and is nothing but a discrete Miura 
transformation [Q. The system ( |3.13| ) with a = —h = i gives an integrable semi- 
discretization of the Gerdjikov-Ivanov equation ( |2.7| ), while ( |3.13| ) with a = h = —1 
gives an integrable lattice version of a higher symmetry of ( p.7[ ) (cf (3.45) in ||27|| ). 
The semi-discrete Gerdjikov-Ivanov system ( |3.13D with b = a* admits the reduction 



= io'q'^*! and not the reduction r° = icrq^_i. 



Through the gauge transformation ( |3.4D with replacing by r°, the linear problem 
]8D and the Lax pair for ( p.l2|) are changed into 



n+l 



and 



zh + 



2 - r'n+lln 



i/2 - 



(3.15a) 



a{z^ - l)h + aq°r°_^ 

1 ^ /yiG ^^G^G^G /1-.G 
"T z'^^n-l'n ^Hn' nHn-l' n-1 


/(;^)(a^g« + fg^i) 


^(-l + ^)(/2-«_i) 

X (azr°_i + f r«) 


K-1 + ^)(^2 - r°g^_i) - ar°g° 
- fe'^^+ig^-i - hr'^^-^q^r^q^_^ 



.(3.15b) 
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Putting the Lax pair ( |3.15| ) into the zero-curvature condition (|2.9|), we exactly obtain 
the semi-discrete Gerdjikov-Ivanov system ( |3.13| ). The L„-matrix ( p.l5a| ) depends on 
and thus it is not ultra-local. It is an open question whether the system 

( p.l3| ) possesses a Lax pair with an ultra-local L„-matrix. In the continuum limit of 
space, we regenerate the eigenvalue problem for the Gerdjikov-Ivanov system. 
By a generalization of the transformation ( |3.14 ) 

for (|2.12| ), we obtain a generalization of (|3.13|) : 

Qn.t - + hqn-i + {a- h)qn + ag„+i(Ar„+i - r„)g„ - &g„(Ar„+i - r„)gn-i 

+ a\qn+irn+iqnrnqn - bXqnrn+iqnrnqn-i = O (3.16a) 
rn,t - hrn+i + ar„_i + (6 - a)r„ - 6r„+i(g„ - Ag„_i)r„ + ar„(g„ - Ag„_i)r„_i 
+ bXrn+iqnrnqn-irn - aAr„g„r„g„_ir„_i = O. 



(3.16b) 

This system is considered as a one-parameter deformation of the matrix Ablowitz- 
Ladik system (|2.12| ). The change of variables g„ = g^^^e'^'^^^^"^'^"'^^'^ ^^^1*, r„ = 
r^A~"e~[*-'*'^^-'"+*^^^'^ ^•'^l* and scalings of a, b cast ( p.l6| ) into ( p. 131) . Transformations 
of this kind are very useful for applications of the discrete DNLS systems and we utilize 
them in sections t4.1|-ff!2. 



3.4- Semi- discrete generalized DNLS system 

In sections pl4|^.5|, we assume that dependent variables are scalar. As a discrete 
analogue of ( p..6|) , we consider the following transformation for the semi-discrete 
Gerdjikov-Ivanov system ( p.l3|) : 

J= — oo J J j= — oo ^ J ^ 

Here 7 is a real parameter. It should be noticed that both log(l — q^-ir^) and 
log(l + q^i^n) conserved densities for ( |3.13|) . The form of the transformation is 
so chosen that the reduction = iag^J^i results in r„ = iaq*_i. 
The inverse of the transformation ( |3.17| ) is written as 

^--^-nrf^ ^»-.n!#^. (3.18) 

9i{-(lj-irj) .^^^ g^[qj_irj_i) 

Here g^^/iy) is defined by the functional equation 

By the successive approximation, gj{y) for a general value of 7 is expressed as the power 
series: 

g^{y) = 1 - 271/ + 7(67 + l)y^ + ■■■. (3.19) 
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Substituting ( |3.18| ) into the semi-discrete Gerdjikov-Ivanov system ( p. 131 ), we obtain an 
integrable lattice version of the generalized DNLS system: 

1 



g-yi-Qrirn+l) 



[1 + (27 + l)qnrn9j{qnrn)] 



+ hqn-ig^{-qn-irn) 



+ qnTn 



[1 - (27 + l)g„r„+i5(^(-g„r.„+i)] 



-g-y{,qnrn) 

+ (a - h)qn - ^-^hqurn+iqug-ii-qnTn+i) - 2'-faqnrnqn9-y{qnrn) = (3.20a) 



rn,t - brn+ig^i-qnTn+l) 



+ avn-ig^iqn-irn-i) 



+ qnTr. 



[1 - (27 + l)qn-irng-y{-qn-irn)] 

1 



- qn-lTn 



1 + (27 + 1) 

qriTng-f 



-9'y[-qn-irnj 

+ (6 - a)r„ + 27ar„g„r„5(^(g„r„) + 276r„g„_ir„5(^(-g„_ir„,) = 0. (3.20b) 

The system ( p.2U| ) with a = —b = i gives an integrable semi-discretization of ( |1.5| ). As 
was expected from the construction of ( |3.17| ), ( p.20| ) with b = a* admits the reduction 
r„ = icrq*_i. Although (|3.20| ) is not an explicit expression for general 7, we can obtain 
approximate expressions by using the power series ( |3.19|) . 
For some choices of 7, g-y{y) is written explicitly, e.g. 



9-i{y) 



l + v/T^ 
^o(y) = 1 



9-i{y) 
9i{y) 



2 



9-iiy) 



1 ^ 1 
1 + 4?/ + 2^ 



1 + v^TT^' 

When 7 = —1/2, (|3.20|) coincides with the semi-discrete Kaup-Newell system (|3.1| ) for 
scalar variables. 



3.5. Alternative semi-discrete Chen-Lee-Liu system 

Setting 7 = —1/4 in ( p. 201) , we obtain an integrable semi- discretization of the Chen- 
Lee-Liu system which is different from (|3.9|) . We write the equation for g„: 

qn,t - a{l + {qnTnf + qnrn\/l + {qnrnY}{[\/l + {qn+lTn+lY + qn+lTn+l] 

X [a/1 + {qnTn+iY - qnrn+i]qn+i - ^n} + b{l + (g„r„+i)^ - g„r„+iv^l + (g„r„+i)2} 
X {[a/1 + (gri_ir„)2 - g„_ir„][A/l + (^nrn)^ + gnr„]g„_i - g„} = 0. (3.21) 

Here we have changed the scaling of g„r„ to eliminate fractions. The alternative semi- 
discrete Chen-Lee-Liu system (( p.21|) and the equation for r„) with b = a* admits 
the reduction = io'q*_i. In this respect, the alternative system is superior to the 
semi-discrete Chen-Lee-Liu system (|0|). 

The system (^) for scalar variables is related to the semi-discrete Gerdjikov-Ivanov 
system ( |3.13| ) via the transformation 



n " 1 

qn=qn n (1 - 'if'^^'j^ ''n=< H TZ^C 



q-_,r- 
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These relations among the semi-discrete DNLS systems ( p.9| ) and ( p.l3| ) can be 

generahzed to the case of matrix- valued variables (cf p9|). 

4. Systems related to the semi-discrete DNLS 

In this section, we obtain integrable semi-discretizations of various systems through 
transformations and reductions for the semi-discrete DNLS systems. 

4-1. Semi- discrete mixed NLS systems 

The so-called mixed NLS equations are obtained as a mixture of an integrable DNLS 
equation and the NLS equation [|1^, |2^, |2^, |2^, |30|, To give an actual example, we 



add +cg^r and — cr^g to the left-hand sides of (|1.5|) , respectively. Since the mixed NLS 
equations are connected with the original DNLS equations through transformations of 
variables, they are also integrable. By a semi-discrete version of such transformations 

for the semi-discrete DNLS systems, we straightforwardly obtain integrable semi- 
discretizations of the mixed NLS systems. 

4-2. Semi- discrete NLS equations and semi- discrete mKdV equations 

For special choices of the parameters 6', yj, the semi-discrete mixed NLS systems are 
reduced to novel lattice versions of the NLS equation and the mKdV equation. 

Substituting = (— l)"e^*^^~°)*g„, = (— l)"e^'^""*)V„ into the semi-discrete 
Kaup-Newell system (|3.1|) , we obtain 

qn,t + - g„+ir„+i)~^g„+i + a(/i - qnTnY^qn - Kh - g„r„+i)"^g„ 

- h{Ii - qn-iTnY^qn-i + 2(6 - a)g„ = O (4.1a) 



rn,t + Kh - rn+iqn)' Tn+i + 6(/2 - r„ - a(/2 - r„g, 

- a(/2 - r„_ig„_i)"V„,„i + 2(a - b)r„, = O. (4.1b) 

Putting q^ = (— l)"e^*^^~"^*g„, = (— l)"e^*^°~'')*r„ into the semi-discrete Chen-Lee-Liu 
system (|3.9| ), we obtain 

qn,t + a{qn+i + qn){h + rnqn) - Kqn + qn^i){h - r„g„_i)"^ + 2(6 - a)g„ = ^^^^ 
rn,t + Kh - r„+ig„)~^(r„+i + r„) - a^h + rnqnWn + Tn-i) + 2(a - 6)r„ = O. 

Substitution of g.^ = (— l)"e^'^^'"^*g„, = (— l)"e^*^"~^)*r„ into the semi-discrete 
Gerdjikov-Ivanov system ( |3.13|) yields 



qn,t + aqn+i - bqn-i + (6 - a)qn + ag„+i(r„+i + r„)g„ - 6g„(r„+i + r„)g„_i 

+ aqn+irn+iqnTnqn - hquTn+iqurnqn^i = O (4.3a) 
rn,t + hrn+i - ar„_i + (a - 6)r„, + 6r„+i(g„ + g„_i)r„ - ar„(g„ + g„_i)r„_i 

+ hrn+iqurnqn-lTn " nqnT nqn-lT n-1 = O. (4.3b) 
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The systems ( [4.1| )-( ^73D give integrable semi-discretizations of the matrix NLS equation 
for a = — 6 = — i and the matrix mKdV equation for a = 6 = 1, respectively. 

In the continuous theory, the matrix NLS equation (|1.1|) and the matrix mKdV 
equation (see (4.24) in ^ or (2.12) in |]^) admit either of the reductions v = au* and 



au\ The system 



with b = a* admits the reduction r„ 



and not 



the reduction r„ = aq _i. This is similar to the property of the matrix Ablowitz- 

" 2 

Ladik system ( |2.12| ) explained in section p.2| . The system ( ^.1| ) with b = a* admits 
the reduction r„ = as well as the reduction r„ = crq*_x. In this respect, the 

differential-difference scheme ( f4.1| ) faithfully inherits the property of the continuous 
hierarchy. This enables us to obtain integrable semi-discretizations of reductions of 
the matrix NLS hierarchy, e.g. the coupled NLS equations H, the coupled Hirota 
equations P, jS^, |5^ and the coupled Sasa-Satsuma equations (cf section 



4-3. Semi- discrete Sasa-Satsuma equations 

If dependent variables are real matrices, the Hermitian conjugation coincides with the 
transposition. The matrix mKdV equation in this case takes the following form 0, Q 



O. 



Ut + Uxxx ~ 3cr(Ux^UU + U^UUx) 

Let us consider the vector reduction: u = (ui, . . . , U2r. 
set of complex variables by ipk = U2k-i + iu2k {k = 1, . 

(cf (4.26) in i) 



Sasa-Satsuma equations [|10|, | 



54, 



'^x-^^mnx^ 



0. 



(4.4) 

e M^'" Defining a 

, m), we obtain the coupled 

(4.5) 



1 and r„ = cr*g„_i gives an integrable semi- 



Here = (^i, V^m) eC™. 

The system ( [4.1|) with a = b 
discretization of the real matrix mKdV equation ([4.4|) . Following the same procedure 
as in the continuous case, we obtain an integrable lattice version of the coupled Sasa- 
Satsuma equations (fl.5|). 



4.4- Semi- discrete KdV equations 

If we consider the reduction a = b = 1, qn = I + Pr, 



I for ([4.1|) , we obtain a matrix 



generalization of a semi-discrete KdV equation |^8|, 

Pn,t 



Pn+1 



Pn-1- 



Through the reduction a = 6 = 1, g„ = — / + 



I for ( [4.3|) , we obtain a matrix 



version of a semi-discrete KdV equation |^ 

yn,t + Vn+lVn " vlVn-X = O. 
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4-5. Semi- discrete Burgers systems 
The matrix Burgers equation, 

- C + = o (4.6) 

is obtained as the reduction r'^ = —I for the Chen-Lee-Liu equation ( p.6|) (after 
scahngs). The Hopf-Cole transformation = —FxF~^ casts the hnear diffusion 



equation Ft — F^x = O into (|4.6| ) 



In a similar way, we obtain integrable semi-discretizations of the Burgers equation 
together with its higher symmetry. Equating with — / in the semi-discrete Chen- 
Lee-Liu system (pl9|), we obtain 

<fn,t - <(ln+l - <fnW ' Qn) ' K<fn ' <fn-lW " = O. (4.7) 

The substitution I — q'^ = Fn+iF~^ relates {^7f\ ) with the linear equation 

- aFn+i + bFn-i + {a- b)Fn = O. 
Setting r„ = — 1 in the alternative semi-discrete Chen-Lee-Liu system ( p.21| ), we obtain 



(ln,t - \/l + q^{a[Jl + - - a[^l + q^ - q^lqn 



+ b[y/l + g2 + _ ^[^17^+ g„-i]g„-i} = 0. (4.8) 
The substitution [a/1 + — g„]^ = fn+ifn^ connects ( |i7^ ) with the linear equation 

fn,t - afn+1 + bfn^l + {a - b)fn = 0. 

The systems ( ^.7| ) and ( ^.8| ) give two lattice versions of the Burgers equation for 
a = —b = 1 and its higher symmetry for a = b = 1, respectively. 

5. Semi-discrete massive Thirring-type models 

In section ^, we obtained the Lax pairs for the semi-discrete DNLS systems. If 
we can replace the M„-matrices with appropriate ones, we obtain integrable semi- 
discretizations of the massive Thirring-type models. For this purpose, the Lax pair 
introduced in [Appendix A| provides useful information. 



5.1. Semi-discrete Kaup-Newell type 



For the semi-discrete Kaup-Newell system in section |3.1| , we replace the M„-matrix 
( p.5b|) with the following one: 



n 1 

Z — - 



lj^_2(^-l)0K^K 2i/(^)^- 



(5.1) 
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Here (j)^ and Xn new variables which are, respectively, li x I2 and I2 x li matrices. 
Substituting the Lax pair ( |3.5a| ) and ( |5.1| ) into the zero-curvature condition ( p.9| ), we 
obtain 



n,t 



An 



o 



An+1 



(5.2) 



i(x^::+x::+i)-2(«x::+x::-,i0X 

o 

ir'' = O. 

The system ( ^.2|) gives an integrable semi-discretization of the massive Thirring model 
of the Kaup-Newell type (see (A. 4) in [^). As we expect from the property of the 
semi-discrete Kaup-Newell system (p.l|), (^.2|) admits either the reduction = iaq^*i, 



K 
An 



i(J(p^*i or the reduction = icrg^Ji , Xn 



1 • 



We can eliminate and to obtain an integrable semi-discretization of the 
Mikhailov model (see (4.21) in [§3] or (3.47) in [g^]), 

€,t - €+i,t + € + €+1 + 2(Cx:^+iC - ^^ix'^a+i^+i) = 
xlt - xl^,,t + xl + xl+, - ^{xlKxl - xl+,<t^lxl^,) = 0. 

Through the reduction = (— l)"(c<j„ + //4), Xn = (~l)"-^5 (|5.3|) collapses into the 
simplest version of the dressing chain 



5.2. Semi- discrete Chen-Lee-Liu type 

For the semi-discrete Chen-Lee-Liu system in section we replace the M„-matrix 
( p. 12 bp with the following one: 

M- = ^\ 2i/(^)0^:? 

" --U7l)(-l + ^)x:^ -l/, + 2(z-i)x^ 

Plunging the Lax pair (|3.12a|) and (|5.4|) into the zero-curvature condition 
obtain 



(5.4) 
9D, we 



ql, + 2i0^ + 2g^^x: 
- 2ix^+i - 2x 



c,,c ,c 



o 



n.t 
AC 



n+l'^n+l'"n 



o 



^n+1 + i^n + <lnr-n€ 



o 



(5.5) 



Xn Xn+1 ^^n Xn+l1: 



o. 



The system ( ^.5|) gives an integrable semi-discretization of the massive Thirring model 
(see (A. 3) in [^). As we expect from the property of the semi-discrete Chen-Lee-Liu 
system ( p.9| ), ( |5.5| ) admits neither the reduction r, 
reduction r!? — --^ — ^-^-'^'^'^ 



An 



lac 



n—k 



nor the 



c t c 

1^9n-fc' Xn 



ic"0n-A;- -'-^ section |5.4| , we obtain an alternative semi- 



discretization of the massive Thirring model without this drawback. 



Integrable discretizations of derivative nonlinear Schrddinger equations 



17 



5.3. Semi- discrete Gerdjikov-Ivanov type 



z — 



^1 



(5.6) 



For the semi-discrete Gerdjikov-Ivanov system in section |3]^, we replace the M„-matrix 
( p.l5b|) with the following one: 

Putting the Lax pair (|3.15a|) and (|5.6|) into the zero-curvature condition ( p.9|) , we obtain 

(lit + i(0n + 0n+l) + 'A(lnrn€ ' ^n+l'^n+l^n) = O 
<,t - i(Xn + Xn+l) + 'A^n^-lXl " Xn+l^O =0 
€ - €+1 + nn + (9n^°0n + = O 

Xn Xn+l ^^n i.^n^n—lXn ~^ Xn+l^n^n^ 

The system ( |5.7| ) gives an integrable semi-discretization of the massive Thirring model of 
the Gerdjikov-Ivanov type (see (A. 4) in with interchanging t and x). As we expect 
from the property of the semi-discrete Gerdjikov-Ivanov system ( |3.13D , ( p.TD admits the 



reduction r° 
If 

model 



If we eliminate '^"'^ '^■■'^ 



i(T</)° *i and not the reduction r° = iaq 



Gf 



1^0! ^ 



and 'we obtain another semi- discretization of the Mikhailov 



- {h + qn+irn+l)(ln,t ' Qn ' Q'. 
^G\^G ^G 



G 

n+1 



o 



(5.8) 



with ([5.3|), we notice the interchange of space and time. Through the 



Comparing 

reduction = (— 1)"(^„ — I), r. 



n 



-1)"J for ( |5.8| ), we obtain the following equation: 



5.4- Alternative semi- discrete Chen-Lee-Liu type 



We assume that dependent variables are scalar and consider the transformation ( |3.17| ) 
together with 

j=—oo ^ J ^ j=—oo ^ J ^ 



Substituting the inverse transformation into (|5.7| ), we obtain an integrable semi- 
discretization of a generalized massive Thirring model. When 7 = —1/2, the obtained 
system coincides with (^.2|) . Setting 7 = —1/4, we obtain an integrable semi- 
discretization of the massive Thirring model which is different from ( ^.5|) . We write 
the equations for g„ and 0„: 



qn,t + i[A/l + {qnTnY + qnTn^n + 'A^ ^ + (g„r„+i)2 - g„r„+i]0n+i 

(5.10a) 

[a/1 + (g„r„)2 + g„rn]0„, - [a/1 + (g„r„+i)2 - g„r„+i]0„+i + ig^ = 0. (5.10b) 
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Here we have changed the scahng of g„r„ to ehminate fractions. Unhke ( |5.5| ), the 
alternative semi-discrete massive Thirring model (( |5.10| ) and the equations for r„ and 
Xn) admits the reduction = iaq *_i , Xn = 



1 • 

'2 '-^2 



6. Concluding remarks 

In this paper, we have investigated integrable discretizations of the DNLS systems from 
the point of view of a Lax-pair formulation. We have found that lattice versions of 
the Kaup-Newell system, the Chen-Lee-Liu system and the Gerdjikov-lvanov system 
are embedded in a generalization of the Ablowitz-Ladik system. With appropriate 
gauge transformations, the eigenvalue problems for the lattice systems coincide with 
the continuous counterparts in the continuum hmit. Using a discrete analogue of the 
phase transformation ( [1.6|) , we obtain a lattice version of the generalized DNLS system 
by Kund. All the discrete DNLS systems but the semi-discrete Chen-Lee-Liu system 
( ^IqD admit the reduction of complex conjugation between two dependent variables. 
This property is indispensable for applications such as difference schemes for numerical 
computation or modelling of nonlinear lattice vibrations. We stress that the discrete 
DNLS systems have relation to a variety of discrete integrable systems. Through changes 
of variables and reductions, we obtain integrable discretizations of the mixed NLS, 
matrix NLS, matrix KdV, matrix mKdV, coupled Sasa-Satsuma, Burgers equations, 
etc. From the eigenvalue problems for the discrete DNLS systems, we derive integrable 
discretizations of the massive Thirring-type models. 

Besides the Ablowitz-Ladik system, there are other integrable discretizations of 
the NLS hierarchy. We conjecture that other discrete DNLS systems are embedded 
in the discrete NLS systems. In fact, we obtain another semi-discretization of the 
Gerdjikov-lvanov system in connection with the semi-discrete NLS system studied in 
|58| (see also |3^, ^). However, the obtained system does not admit the reduction of 



complex conjugation. Thus it is less attractive than ( p.l3| ). The example implies that 
the Ablowitz-Ladik formulation is an optimal starting point for the study of discrete 
DNLS systems. 

The NLS and DNLS systems for scalar variables possess tri-Hamiltonian structure 
|59| , |60[| . The three sets of Poisson brackets are, respectively, the 5\x — y) type for 
the Kaup-Newell system, the 5{x — y) type for the Chen-Lee-Liu system and the 
6{x — y) type for the NLS system. We should note that all these systems are connected 
via transformations of the dependent variables. We infer that the lattice versions of 
these systems studied in this paper possess tri-Hamiltonian structure. The discrete 
counterpart of the first Hamiltonian structure is ( |3.6| ) for the semi-discrete Kaup-Newell 
system ( |3.1| ). The Ablowitz-Ladik system ( |2.12|) has a set of ultra- local Poisson brackets 



I|, |63[ which is the discrete analogue of the third Hamiltonian structure. Since 
the lattice systems are each connected with the others through changes of variables 
(see section |]), they are at least bi-Hamiltonian. We have not found the discrete 
counterpart of the second Hamiltonian structure for either the semi-discrete Chen-Lee- 
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Liu system ( p.9| ) or the alternative system ( p.21| ). It remains an unsolved problem to 
prove explicitly that the Ablowitz-Ladik system and the semi-discrete DNLS systems 
are tri-Hamiltonian. 
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Appendix A. Semi-discrete sine-Gordon equation as a negative flow of the 
Ablowitz— Ladik hierarchy 

In this appendix, we prove that a semi-discrete sine-Gordon equation and its 
generalization are integrable via the ISM based on the Ablowitz-Ladik eigenvalue 
problem. We introduce the following form of the Lax pair: 



Z Un 
1 



Cn. ZQjYi 



(A.l) 



Here z is the spectral parameter. u„, f„, a„, 6„, c„ are scalar variables. The M„-matrix 
is simpler than the one by Ablowitz and Ladik []TB|, while the L„-matrix is the same. 
Substituting ( [A.l| ) into the zero- curvature condition ( |2.9| ), we obtain the following six 
equations: 

an - an+l + UnCn " t^n&n+l = (A.2a) 
Un{Cn - Cn+l) + Vn{hn " &n+l) = (A. 2b) 

6„ - - M„(a„ + a„+i) = (A.2c) 
c„-c„+i + t;„(a„ + a„+i) =0 (A. 2d) 

Un,t + bn- Unttn = (A.2e) 
Vn,t - Cn- Vnan = 0. (A.2f) 

Thanks to ( [A.2c| ) and ( [A.2d|) , we can write u„ and f„ in terms of a„, 6„, c„: 

bn ~ bn+l Cn — Cn+l 



Ur. 



Then (|A.2b|) is automatically satisfied. Putting ( |A.3|) into (|A.2a|) , we obtain 

al + bnCn = Tj^ 

where 77 is a constant. 
If we set 

an = r] cos On bn = Cn = r] sin 6'„ 
Un and Vn are expressed as 



(A.3) 



(A.4) 



Ur. 



sin 9n - sin 9n+l . (On- 9n+l 
= tan 

cos On + COS On+1 V 2 
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Then ( |A.2^ ) and ( [A.21| ) are combined into a semi-discrete sine-Gordon equation |6^, |6 

On+u - dn,t = ?7(sin6'„ + sin6'„+i). (A. 5) 

By defining a„ by 6'„ = (a„ + ari+i)/2, we rewrite ( |A.5| ) as |2^ 

an+i,t - 0!n,t = 2?7 sm j (A. 6) 

under appropriate boundary conditions. 

More generally, 6„ and c„ are functionally independent. Substituting a„ = 



\/ 'ff ~ bnCn into ( |A.3|) , ( |A.2e| ) and ( |A.21|) , we obtain a semi-discretization of the reduced 
equation for the 0(4) nonlinear cr-model 



K\/lf - bn+lCn+l + bn+l\/rf - bnC„ 



(A. 7a) 



I H , = , ^ = 0. (A. 7b) 

a///^ - 6„c„ + A/r?2 - bn+lCn+l J ^ - &nC„ + y^r]^ - bn+lCn+1 

If we set rj = and 6„ = — c„ = e*^", ( |A.7|) is simplified to a semi-discrete Liouville 
equation HBBj , 



We have shown that the semi-discrete sine-Gordon equation ( |A.5| ) and its 
generalization (|A.7| ) are associated with the Ablowitz-Ladik eigenvalue problem. Thus 
they are integrated by the ISM straightforwardly. The semi-discrete sine-Gordon 
equation ([A.5|) is associated with another ISM-solvable eigenvalue problem |65 



It is now clear that the difference between the two eigenvalue problems is not essential 



(see also [0]). We obtain a lot of knowledge on the semi-discrete sine-Gordon equation 
and its generalization from the study of the Ablowitz-Ladik hierarchy. 

In section |], we investigate integrable semi-discretizations of the massive Thirring- 
type models. The M„-matrix in ( |A.1| ) with the constraint ( [A.4| ) gives helpful information 
for that. 



Appendix B. Time discretizations of semi-discrete sine-Gordon equation 

In this appendix, we discuss time discretizations of the semi-discrete sine-Gordon 



equation and its generalization obtained in Appendix A . We consider a time 
discretization of (|2.8|) : 

(B.l) 



n+l 



L \& 



^ r 



V 



where the tilde denotes the step-up shift (/—*>/ -|- 1) in the discrete time / G Z. The 



compatibility of (|B.1|) leads to a full-discrete version of the zero-curvature condition 



(B.2) 
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As a discrete-time analogue of ( |A.1| ), we introduce the following form of the Lax pair: 



z 



Ur, 



1 
1 



+ 



h 



(B.3) 



Here h denotes the difference interval of time. Substituting ( B.3 ) into the zero-curvature 
condition ( |B.2| ), we obtain the following six equations: 

(B.4a) 
(B.4b) 
(B.4c) 
(B.4d) 
(B.4e) 
(B.4f) 



"UnCn ~ + Vnbn — V-ribn+l = 

bn ~ bn+1 — Und-n ~ UnCtn+1 = 
Cra — C„_|_i + VnO'ii + VnO-n+l = 
Un-Un + h{bn - Unan) = 
Vn-Vn- h{Cn + ^nCtn) = 0. 

Thanks to ( p.4cj ) and (|B.4e|) , we can write Un and in terms of a„ and 6„: 

bn — bn+l + hanbn+l _ bn — bn+1 ~ httn+lbn 



Ur 



O'n + o^n+i ~ ha^an-^i 



Ur. 



Q-n + 0-71+1 ^ ha^an-^i 



Similarly, using (|B.4d|) and (p.41]) , we express Vn and Vn in terms of a„ and c„ 



(B.5a) 



(B.5b) 



Oft + O-n+l ~ hanttn+l Q-n + O-n+l ~ ha^an+l 

It is easily seen that (p.4b|) is automatically satisfied. Putting (|B.5D into (p.4a| ), we 
^''1 — han) or, equivalent ly 



obtain + 6„c„ = 77 
(a„ + 



^)' + &nCn = r72[l + ^ — 



(B.6) 



Here is a constant. 

The relations (p.5|) lead to the following system. 



bn — bn+l + hanbn+l _ bn — bn+1 — hUn+lbn 
O-n + Cin+1 ~ hdnCln+l Cln + ^n+l ~ hO-n^'n+l 

Cn ~ Cn+1 + httnCn+l Cn ~ Cn+1 ~ httn+lCn 
O'n + (in+1 ~ handn+1 Cln + Cln+l ~ han^n+l 

where a„, 6„, c„ are related through ( |B.6| ). Substituting 



(B.7a) 
(B.7b) 



2 +r/Wl + (y)^cos^„ 



c„ = ?7\/l + (^)^sin6'^ 



into (p.7|) , we obtain an integrable time discretization of the semi-discrete sine-Gordon 
equation (|A.5|) . More generally, putting 



+ 



b c 



2 V 2 ' 

into (|B.7|) , we obtain an integrable time discretization of the semi-discrete reduced 
nonlinear cr-model ( [A.7|) . 
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We shall extract a celebrated full-discrete sine-Gordon equation from ( |B.4| ). The 
equation is interpreted as an integrable time discretization of ( [A .61 ). We introduce the 
following parametrization: 

^V"^ , „ /7~7^a2„„„ + (B8a) 

(B.8b) 

(B.Sc) 



+ ^Y1 + (y) cos(^- 



bn = Cn = + {-^f ^^'^( 



tan 



a. 



hri-,2 ■ fan + an+i 
smi 

an+2 



It is sufficient to consider ( |B.4aD , ( |B.4cD , ( |B.4e[ ) among ( |B.4| ). For brevity, we employ 
the following abbreviations: 



an + 



2 +'"'Vi + 't'^- 



4 4 
Under the parametrization (|B.8| ), we rewrite (p.4a| ), ( |B.4c| ), (|B.4e|) respectively as 

cos Xn COS Un+l COS Un COS 



(B.9a) 
(B.9b) 
(B.9c) 



cos(x„ - yn+i) cos{y 

k sin Xn cos yn+i + cos Xn sin yn+i sin yn cos -|- k cos sin Xn+i 
cos(a;„ - y„+i) cos(?/„ - 

A; sin Xn cos — A;~^ cos Xn sin j/^+i sin ?/„ cos x„+i — cos ?/„ sin 
cos(x„ - yn+i) cos{y 

With the help of ( [B.9aD , we combine ( |B.9b| ) and ( [B.9c| ) into one equation: 

ktanXn = tanyn- (B.IO) 

Conversely, all of the relations ( p.9|) hold if x„ and yn satisfy (P.10| ). Equation ( p.lCI|) is 
the celebrated full-discrete sine-Gordon equation |6^. It is easily cast into the standard 

'On + On+l + ftn + "n+l" 



form 19 



sm 



an + an+i - ar. 



hri 
~2 



: sm 



We have obtained integrable full-discretizations of the sine-Gordon equation and 
the reduced nonlinear cr-model. The obtained equations are integrable via the ISM 
based on the Ablowitz-Ladik eigenvalue problem. 
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